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ABSTRACT 

The  following  is  a  report  on  the  design  and  construction  of  a 
fixed-base  variable-stability  simulator  facility  combined  with  the  task 
of  landing  to  a  carrier  or  runway. 

The  solution  was  mechanized  on  a  hybrid  computer  with  the  analog 
computer  solving  the  equations  of  motion  and  the  digital  computer  used 
for  storage,  control  and  graphics  generation.  The  display  was  in  the 
form  of  a  computer-drawn  picture  on  a  graphics  terminal.  Control  was  by 
a  simulated  cockpit  placed  in  front  of  the  display  and  connected  to  the 
analog  computer. 

Dynamic  validation  was  considered  excellent  with  the  modal  periods 
of  the  simulated  aircraft  agreeing  very   closely  with  those  of  the  actual 
aircraft. 

The  visual  display  was  deemed  very  good  as  sufficient  visual  cues 
were  provided  to  enable  consistent  landings  by  experienced  pilots. 

This  project  was  undertaken  not  as  a  design  of  training  aid  but 
rather  as  a  research  tool  for  further  studies  in  control  systems,  human 
engineering  and  aircraft  dynamics. 
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I.   INTRODUCTION 

The  purpose  of  this  project  was  to  construct  a  multiple  task  simu- 
lator, an  inherent  quality  of  which  had  to  be  the  ease  of  its  extension 
as  a  research  tool  to  various  real  world  problems  that  arise.  In  addi- 
tion, the  use  of  interactive  computer  graphics  would  be  studied  as  a 
means  of  visual  display  for  the  simulator.  The  decision  was  made  to 
design  the  simulator  to  the  initial  task  of  approach  and  landing  to  an 
aircraft  carrier  deck. 

The  hub  of  the  simulator  system  is  a  hybrid  computer,  part  digital, 
part  analog.  The  analog  computer  solved  the  equations  of  motion  while 
the  digital  computer  was  used  for  control,  storage  and  graphic  build-up. 
The  display  was  handled  by  a  graphics  processor  and  terminal  connected 
directly  to  the  digital  computer.  Control  of  the  simulator  was  through 
a  simulated  cockpit  placed  in  front  of  the  graphic  display,  with  stick 
and  throttle  outputs  tied  directly  to  the  analog.  This  then  formed  the 
control  loop  as  diagrammed  in  Figure  1. 

At  the  beginning,  a  decision  on  the  allocation  of  equipment  func- 
tions had  to  be  made.  Several  all -digital  moving-base  simulators  were 
already  constructed  [Ref.  1].  As  space,  money  and  time  were  limited, 
moving  base  was  out  of  the  question  but  a  pure  digital  solution  was  not. 
On  the  other  hand  the  department  already  had  a  "Mini  Link"  Simulator 
utilizing  an  all  analog  solution.  However,  due  to  the  lack  of  non-linear 
equipment  and  the  need  for  generating  a  visual  perspective  display,  the 
pure  analog  solution  could  not  be  implemented.  Then  too,  the  digital 
solution  was  too  slow  for  the  needed  integration  and  did  not  provide 
for  a  pilot-station  computer  interface.  The  obvious  solution  was  the 
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Figure  1.  Lines  of  Communications 
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hybrid  computer.  The  non-linearities  of  the  aircraft  coefficients,  the 
resolution  of  the  forces  and  moments  and  the  generation  of  the  display 
would  be  tasked  to  the  digital  computer.  The  integration  of  the  equa- 
tions of  motion  and  the  interfacing  would  be  done  by  the  analog.   In  this 
manner  the  best  features  of  both  computers  could  be  utilized. 

The  display  would  be  a  computer-drawn  picture  on  a  graphics  terminal 
connected  to  the  digital  computer.  Construction  of  the  picture  would  be 
done  by  the  digital  computer  while  the  display  and  refresh  would  be  done 
by  the  graphic  digital  processor. 

The  "cockpit"  was  simply  a  seat  and  attached  controls  placed  in 
front  of  the  graphics  CRT.  Stick  and  throttle  outputs  from  the  cockpit 
were  then  fed  directly  into  the  analog  computer,  thus  closing  the 
control  loop. 

The  design  features  incorporated  into  the  Carrier  Approach  Landing 
Simulator  (CALS)  were: 

1.  Ability  to  change  the  aircraft  being  simulated 

2.  Capability  of  simulating  non-linear  aerodynamic  data 

3.  Automatic  scaling  of  the  analog  computer  solution 

4.  Fixed-base  inside-out  display 

5.  Ability  to  expand  or  change  the  simulator  task 

6.  Simplicity  of  setup  and  operation 
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II.  THE  SIMULATOR 

The  simulator  is  capable  of  accurately  simulating  vehicle  dynamic 
response  in  conjunction  with  a  dynamic  inside-out  display  with  its  view- 
ing plane  affixed  to  the  vehicle.  The  analog  solution  is  updated  at  a 
rate  of  18  samples  per  second  which  leads  to  excellent  dynamic  response. 
A  new  graphic  display  image  is  generated  at  the  rate  of  18  frames  per 
second  with  a  digital  processor  refresh  rate  of  40  frames  per  second. 
This  leads  to  a  display  with  virtually  no  flicker  and  in  general  yery 
good  graphic  qualities. 

The  simulator  was  designed  to  accept  the  following  inputs  in  the 
form  of  punched  cards:  the  coordinates  of  the  runway;  the  runway  cons- 
tants (glide  slope  and  lens  placement);  the  aircraft  constants  (eight, 
moment  of  inertia,  wing  area  etc.);  initial  conditions  of  the  problem 
(initial  position,  altitude  and  velocity)  the  maximum  expected  range 
of  the  state  and  control  variables  in  the  form  of  scale  factors:  and  the 
aerodynamic  coefficients  of  the  simulated  aircraft.  The  preparation  of 
this  data  deck  is  outlined  in  Appendix  F. 

The  simulator  then  processes  this  data,  autoscales  the  analog 
computer  based  on  the  given  and  self-generated  scale  factors,  sets  the 
potentiometers  of  the  analog  and  places  the  initial  conditions  on  the 
proper  integrators.  Control  of  the  simulator  is  then  transferred  to 
the  pilot  who  responds  as  per  the  displayed  instructions.  At  the 
completion  of  a  run  the  simulator  analyzes  the  landing  and  displays  the 
results.  This  analysis  includes  sink  rate  at  touch  down,  line-up, 
airspeed  and  type  of  landing  (crash,  bolter  or  arrested  landing  on  a 
certain  wire).  At  this  time  options  are  offered  to  either  fly  again 
with  the  same  parameters,  change  the  parameters  or  stop  the  program. 
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Besides  these  options,  the  operator  can  suppress  printed  output  if 
desired,  select  a  yaw  damper  (ON/OFF)  and  inhibit  translational  movement 
while  flying  the  aircraft.  Also,  at  the  disposal  of  the  operator,  is  an 
eight-channel  strip  recorder  and  an  oscilloscope  that  are  patched  into 
the  analog  computer.  This  gives  the  operator  the  ability  to  output  the 
time  history  of  any  variable  he  wishes. 

The  cockpit  consists  of  a  chair  upon  which  the  throttle  and  control 
stick  have  been  mounted.  The  stick,  mounted  on  the  right  arm,  is  a 
Gemini  control  stick  providing  all  three  controls:  yaw,  pitch  and  roll. 
At  the  top  of  the  stick  is  located  one  of  the  control  buttons  used  in 
the  control  of  the  program  from  the  pilot  station.  Mounted  on  the  left 
arm  is  a  throttle  plate  containing  a  small  throttle  and  another  control 
button.  The  outputs  of  the  cockpit  are  fed  via  a  wire  bundle  to  the 
rear  of  the  analog  computer  for  processing. 

A  complete  operating  manual  is  contained  in  Appendix  H. 
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III.  SIX  DEGREE  OF  FREEDOM  AIRFRAME  EQUATIONS 

R.M.  Howe  [Refs.  2  and  3]  discusses  three  coordinate  systems  -  body, 
wind,  and  stability  axes  -  and  concludes  that  the  use  of  wind  axes  for  the 
translational  equations  and  of  body  axes  for  the  rotational  equations  of 
motion  is  the  best  choice.  As  can  be  seen  from  Figure  2,  the  body  axes 
differ  from  the  stability  axes  by  the  aircraft  angle  of  attack,  while 
the  stability  axes  differ  from  the  wind  or  flight  path  axes  by  the  side 
slip  angle. 

Assuming  that  the  three  angular  rates  of  the  aircraft  (P,Q,R) 
about  the  body  axis  are  given,  the  angular  rates  about  the  stability 
axis  (P  ,Q  ,R  )  can  be  computed  (Eqn  101-103). 

P$  =  P  cos  a  +  R  sin  a  (101) 

Qs  =  Q  (102) 

R$  =  P  sin  a  +  Q  cos  a  (103) 

These  angular  rates  can  then  be  normalized  using  equations  104-106. 
P  b 

ps  ■  w-  (104) 

Q  b 
Qs  =  gf  (105) 

R  b 
Rs  =  fr  (106) 

The  aerodynamic  forces  expressed  along  the  stability  axes  can  then 
be  calculated  (Eqn  107-109). 

Fx      /-scn\ 

aero  m       i   _D  (107) 


m     n   \    m 
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Figure  2.  Aircraft  Coordinate  Axis 
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The  forces  in  the  body  axes  are  then  resolved  into  forces  in  the 
stability  axes  (Eqn  110-112).  Thrust  and  the  force  of  gravity  are 
included  at  this  time.  Thrust  was  assumed  to  have  only  an  X  component 
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The  forces  are  then  resolved  from  the  stabil ity  axes  to  the  wind  or 
flight  path  axes  (Eqn  113-115). 
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From  these  equations  the  translational  equations  of  motion  can  be 
determined.  Solving  for  the  derivatives  yields  equations  116-118. 
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-3  = 
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(118) 


Further,  the  Euler  angular  rates  can  then  be  calculated  (Eqn  1 19- 


121). 
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Now  focusing  our  attention  on  the  moment  equations,  the  three  equa- 
tions about  the  stability  axes  can  be  written  as  follows  (Eqn  122-124). 
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These  moments  are  then  resolved  into  the  body  axes  yielding 

Equations  125-127. 
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The  rotational  equations  of  motion  can  then  be  written  (Equations 
128-130). 
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As  suggested  by  Howe  [Ref.  2],  a  simplified  method  for  computing 
the  velocities  in  the  inertial  frame  yields  equations  131-133.  This 
approach  uses  small  angle  approximations  in  its  development. 

S  =  V  cos  e  cos  $  (131) 

As 

S  =  V  cos  e  sin  i>  (132) 

5  =  V(-sin  e  +  sin  a  cos  8  cos  <j>)  (133) 

The  various  derivatives  that  were  solved  for  in  the  above  equations 
can  then  be  integrated  yielding  the  state  variables  that  were  assumed  to 
exist  at  the  start  (Eqns  134-145) 

P  =  /P  dt  (134) 

Q  =  /Q  dt  (135) 

R  =  /R  dt  (136) 

a  »  /a  dt  (137) 

6  =  /3  dt  (138) 
<p  =  fi  dt  (139) 
e  =  si  dt  (140) 
<f>  -  y<t»  dt  (141) 
V  =  /V  dt  (142) 
x  =  /Sv  dt  (143) 
y  =  /S  dt  (144) 
z  =  /Sz  dt  (145) 
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The  equations  themselves,  when  examined,  seemingly  divide  themselves 
into  two  groups.  One  that  can  be  solved  by  digital  means  and  one  that 
can  be  solved  by  the  use  of  the  analog  computer.  In  fact  this  logical 
division  is  the  one  chosen  for  the  simulator. 
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IV.  GRAPHICS  PRESENTATION 

Since  the  primary  scan  during  a  carrier  approach  consists  of  only 
three  things,  the  "meatball"  or  the  Fresnel  Lens  Optical  Landing  System, 
the  aircraft's  angle  of  attack  and  the  line-up  with  the  runway  center- 
line,  it  was  decided  to  incorporate  only  these  three  visual  cues  into 
the  display. 

The  angle  of  attack  was  represented  by  an  angle  of  attack  indexer 
situated  in  the  upper  right  hand  corner  of  the  display.   Its  operation 
is  exactly  the  same  as  the  one  found  in  the  jets/turboprops  operating  in 
the  Navy  today.  The  indexer  has  five  descrete  states  as  shown  in 
Figure  3. 

v  v 

D  □  □ 

A   A 

States     12       3        4        5 

Figure  3.  Angle  of  Attack  Indexer  States. 
Optimum/Landing  angle  of  attack  was  represented  by  State  3.  As  one 
increases  the  angle  of  attack  of  the  aircraft,  State  2  would  appear 
indicating  the  aircraft  was  a  "little  slow"  and  eventually  State  1  would 
appear  indicating  that  the  aircraft  was  "slow".  Proceeding  the  other 
direction  from  State  3,  that  is  decreasing  the  angle  of  attack,  first 
State  4  would  appear  indicating  the  aircraft  was  a  "little  fast"  and 
finally  State  5  would  appear  indicating  the  aircraft  was  "fast".   In 
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terms  of  degrees,  assuming  optimum  angle  of  attack  to  be  0°,  the 
following  would  exist: 

STATES  RANGE  OF  AOA 

1  Greater  than  2.0° 

2  Between  1 .0°  and  2.0° 

3  1.0°  to  -1.0° 

4  Between  -1 .0°  and  -2.0° 

5  Less  than  -2.0° 

The  Fresnel  Lens  is  represented  by  two  fixed  rectangular  boxes  for 
the  datum  lights  and  a  movable  square  for  the  "meatball"  as  shown  in 
Figure  4. 


□ 


Figure  4.  Centered  Meatball. 
The  operation  is  again  the  same  as  found  in  the  Navy  systems.   If 
the  "meatball  is  above  the  datum  bars,  as  in  Figure  5,  the  aircraft  is 
above  the  glide  slope. 


□ 


Figure  5.  High  Meatball . 
If  the  "meatball"  is  below  the  datum  bars,  as  in  Figure  6,  the 
aircraft  is  below  the  glide  slope. 


D 


Figure  6.  Low  Meatball 
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The  object  is  to  keep  the  "meatball"  centered  within  the  datums. 
The  "meatball"  has  a  1.5  degree  visibility  arc  centered  on  the  optimum 
glide  slope  and  will  disappear  if  this  is  exceeded.  Also,  once  you  fly 
past  the  lens  itself  the  "meatball"  will  disappear.  However,  in  this 
simulation  the  datums  will  remain  visible  at  all  times,  only  the 
"meatball"  disappears. 

Since  the  lens  is  affixed  to  the  carrier/runway  in  real  life,  the 
size  of  the  presentation  will  change.  However,  since  the  vast  majority 
of  flight  cues  are  absent  in  the  simulator,  the  size  of  the  simulated 
lens  will  remain  fixed.  This  presentation  is  placed  in  the  upper  left 
hand  corner  of  the  display. 

The  line-up  cue  is  received  from  the  picture  of  the  runway  itself. 
The  picture  shown  is  as  if  a  movie  camera  were  affixed  to  the  nose 
(X-axis)  of  the  aircraft.  The  field  of  view  is  a  square  with  an 
angular  limit  of  plus  or  minus  18.5  degrees  in  the  vertical  and 
horizontal  directions.  Within  this  window  lies  the  horizon  and  the 
runway  outline.  Tic  marks  at  0°  and  plus  and  minus  10°  of  pitch  are 
given  along  the  left  side  of  the  window.  As  the  aircraft  maneuvers 
through  space  the  runway  and  horizon  move  dynamically  as  if  filmed 
by  the  movie  camera.  The  complete  display  is  shown  in  Figure  7. 
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V 

□ 

A 


Figure  7.  The  Complete  Graphics  Display, 
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V.  PROCESSING  THE  DISPLAY 

The  processing  of  the  angle  of  attack  indexer  is  straightforward  as 
it  is  solely  a  function  of  the  angle  of  attack  of  the  aircraft.  A  simple 
test  of  the  angle  of  attack  will  yield  one  of  the  five  previously 
mentioned  states. 

The  Fresnel  Lens  or  "meatball"  is  again  a  simple  task  to  process. 
Knowing  the  x-coordinate  and  the  z-coordinate  of  the  aircraft  and  the 
x-coordinate  of  the  lens,  the  glide  slope  angle  can  be  calculated. 


(xL,zL) 

Figure  8.  Glide  Slope  Computation. 
The  angle  GS  in  Figure  8  is  the  glide  slope  angle  while  (X  ,Z  )  and 
(X.  ,Z. )  are  the  coordinates  of  the  plane  and  lens  respectively. 

GS  =  tan"1((Zp-ZL)/(Xp-XL)) 

Angle  GS  is  then  calculated  and  compared  to  the  optimum  glide  slope, 
and  the  "meatball"  is  drawn  accordingly.   If  GS  is  greater  than  the 
allowable  limits,  or  if  X  <  X, ,  the  "meatball"  is  not  displayed. 

The  processing  of  the  runway  is  by  far  the  most  complicated.  How 
to  transform  a  three-dimensional  object  like  the  runway  into  a  two- 
dimensional  object  capable  of  being  displayed  on  the  screen  of  a  CRT? 
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The  answer  was  found  in  two  works,  one  by  L.G.  Roberts  [Ref.  3]  and  the 
other  by  R.B.  Desens  [Ref.  4].  Roberts  introduces  the  idea  of  a  homo- 
geneous coordinate  system  which  introduces  a  fourth  coordinate  or 
variable  scale  factor.  He  also  constructs  a  single  transformation 
matrix  for  all  points.  This  transformation  matrix  (H-Matrix)  is  a 
product  of  five  matrices:  rotation,  translation,  perspective,  another 
translation  and  scale.  Each  point  is  post-multiplied  by  the  H-Matrix 
and  then  the  scale  factor  is  divided  out  yielding  the  display  coordinates 
Desens  used  this  to  process  pictures  of  an  aircraft  carrier  from  various 
positions.  He  also  reduced  the  H-Matrix  to  a  product  of  three  individual 
matrices  as  shown  in  Figure  9. 


H  = 


M-i    Mq    /»-j    U 

B-i  B^  Bq  0 
XT  YT  ZT  1 


Object  Axis 
Rotation 
(Euler  Angles) 
Viewing  Plane 
Translation 


AA-,  AA2  AA3  0 

BB,  BB2  BB3  0 

CC-,  cc2  cc3  0 

0  0    0  1 


Viewing  Plane 
Orientation 


1  -Y  /F  -Z  /F  -S/F 


0  1 
0  0 
0   Y_ 


Perspective 
Scale 
Offset 


Figure  9.  H-Matrix 
Figure  10  has  the  two  coordinate  systems,  that  of  the  earth  and 
that  of  the  aircraft.  The  viewing  plane  or  camera  is  placed  on  the 
nose  of  the  aircraft  while  the  runway  is  oriented  in  the  object  or 
earth  reference  system.  The  first  matrix  is  made  up  in  part  by  the 
direction  cosines  of  the  Euler  angle  rotation  of  the  object  (runway) 
In  this  case  the  runway  is  fixed.  The  Euler  angle  direction  cosines 
are  given  below  and  are  based  on  the  rotations  shown  in  Figure  11. 
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Figure  10.  Graphic  Coordinate  System, 

Euler  Angle  Direction  Cosines 
A-,  =  cos  B  cos  a 
Ap  =  cos  3  sin  a 
A3  =  -sin  3 

B-,  =  -cos  y  sin  a  +sin  y   sin  6  cos  a 
B?  =  cos  a  cos  y  +sin  y  sin  6  sin  a 
B3  =  sin  y  cos  6 

C-.  =  sin  y  sin  a  +cos  y  sin  6  cos  a 
Co  =-sin  y  cos  a  +cos  y  sin  s  sin  a 
Co  =  cos  B  cos  y 


Order  of  Rotation 
a,  B,  Y 


3* 


G 


Figure  11.  Euler  Angle  Rotation 
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Xj,  Y-j.  and  Z-j-  are  the  translation  distances  between  the  two  axis 
systems. 

X  =  X      -  X 
T    object    view  plane 

Y  =  Y      -  Y 
T    object    view  plane 

T    object    view  plane 

The  second  matrix  is  made  up  entirely  of  the  viewing-plane-orienta- 
tion  direction  cosines.  The  direction  cosines  are  qiven  below  in  terms 
of  the  Euler  angles  of  the  aircraft. 

AA,  =  cos  4>   cos  e 

AAp  =  sin  \i>  cos  <\>  -cos  $  sin  e  sin  <}> 

AA3  =  sin  \p  sin  <j>  +cos  ty  sin  e  cos  <}> 

BB,  =  -cos  e  sin  \\> 

BBp  =  cos  ty  cos  <j>  +sin  e  sin  \\>  sin  <\> 

BB^  =  cos  i>  sin  <j>  -sin  4>  sin  0  cos  <f> 

CC-,  =  -sin  0 

CCp  =  -cos  e  sin  <£  Order  of  Rotation 

CC3  =  cos  e  cos  $  i>   (yaw),  e  (pitch),  <\>   (roll) 

The  third  matrix  consists  of  the  offset  option,  Z  and  Y  ,  which 
allows  one  to  look  at  an  offset  portion  of  the  viewing  plane.  This 
option  is  not  used  on  this  display.  The  matrix  also  contains  the  scale 
factor  (S)  and  the  focal  length  (F).  The  scale  factor  is  set  at  1/2 
and  the  focal  length,  the  distance  between  the  viewer  and  the  viewing 
plane,  is  set  at  1.5  ft.  The  field  of  view  is  controlled  by  S  and  F 
and  is  equal  to  tan"  (S/F). 

A  point  on  the  object  (X,Y,Z)  is  converted  to  homogeneous  coordinates 
by  the  addition  of  another  scale  factor  (W)  yielding  (X1,  Y',  Z',W) 
where  X'  =  WX,  Y'  =  WY,  V   =  WZ.  This  then  is  post  multiplied  by  the 
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H-Matrix  yielding  (X",  Y",  Z",  W).  The  display  coordinates  of  the 
processed  point  can  then  be  found  by  dividing  by  W 

V"  =  Y"/W 

Z'"  =  Z'VW 

The  coordinate  X'"  =  X '  '/W  is  an  indicator  of  depth  and  can  be 
used  for  depth  cues,  however,  it  is  not  used  in  this  display.  Each 
point  of  the  runway  is  passed  through  the  H-Matrix  yielding  the  display 
coordinates  of  each  point  in  turn. 

Turning  to  the  problem  of  the  horizon  and  referencing  Figure  12 
below  the  following  computations  are  presented. 


Figure  12.  Horizon  and  Earth. 
For  an  aircraft  at  point  A  with  altitude  z  above  the  earth  with 
radius  r  ,  the  true  horizon  is  at  point  B.  This  is  depressed  a  distance 
z  from  the  plane  DE  which  is  tangent  to  the  earth  at  point  D.  However, 
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a  pseudo  horizon  exists  at  point  E  which  always  has  a  z  coordinate  of 
zero  and  visually  coincides  with  the  true  horizon  but  has  the  advantage 
of  being  quicker  to  draw. 

From  Figure  12,  since  triangle  ADE  and  ABC  are  similar  this  yields 

DE  =  a_  ..  BC   ^e 
AD  ""  z   AB   AB 

Re(z) 

a  =  ~AB~ 

(z  +  re)2  =  AB2  +  re2 

2        ] 
AB  =  {zd   +  2(re)z)2 

re(z)2 
a  =  =- 


1 

(z  -  2r  )7 
e' 


assuming  2r  to  be  wery   much  greater  than  z  yields: 

1  ! 

a  =  - r  =  K(z)2   where  K  =  (41'8)2^  x  103 

(2re)2 

The  depression  angle  of  the  horizon  is: 

-1  1    1 

dep  =  90°  -  tan  '(a/z)  =  90°  -  tan" ' (K(z)2/z) 

Referring  to  Figure  13,  which  is  a  view  of  the  window  and  horizon, 
the  object  is  to  calculate  YMID,  the  y  intercept  of  the  horizon.  The 
distance  PITCH  can  be  calculated  knowing  the  pitch  angle  and  the  field 
view. 

PITCH  =  e/Field  of  View 
TPITCH  is  PITCH  minus  the  depression  angle  DEPA.  Simple  trig  yields 
cos  <f>  =  TPITCH/YMID,  where  <|>  is  the  roll  angle  of  the  aircraft. 
YMID  =  TPITCH/cos  * 
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Knowing  YMID  and  the  roll  angle  the  horizon  can  then  be  constructed 
Upon  the  completion  of  the  runway  and  the  horizon,  the  lines  are  placed 
through  a  software  window  algorithm  that  will  cause  only  the  lines  with- 
in the  window/field  of  view  to  be  displayed.  This  computation  is  done 
for  every   iteration  of  the  problem  solution. 


Figure  13.  Calculating  Horizon  Intercept, 
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VI.  THE  PROGRAM 

As  stated  in  Section  III,  a  logical  division  of  tasks  exist.  To 
integrate  the  derivatives  of  the  twelve  state  variables  by  digital 
algorithms  takes  too  much  computational  time,  therefore  the  logical 
thing  to  do  was  to  D/A  these  variables  and  integrate  them  on  a  real- 
time basis  using  the  analog  computer.  The  results  are  then  A/D'd  back 
to  the  digital  computer  for  calculation  of  the  derivatives.  This  forms 
the  computational  loop. 

When  a  real-time  simulation  is  undertaken  the  time  involved  for  the 
solution  of  the  problem  must  be  a  minimum.  Great  pains  were  taken  to 
insure  that  this  was  indeed  done.  The  following  is  a  tabulation  of  the 
significant  time-saving  steps  that  were  involved: 

1.  Resolution  of  the  angles  in  the  digital  computer  using  a 
sine  table  lookup. 

2.  Fast  linear  interpolation  used  in  both  the  sine  table 
lookup  and  the  aerodynamic  coefficient  lookup. 

3.  Preprocessing  of  the  aerodynamic  data  to  yield  usable 
aerodynamic  coefficients. 

4.  General  and  liberal  use  of  program  constants  to  combine 
two  or  more  arithmetic  operations. 

5.  Use  of  the  SQRT  function  vice  the  exponential  function. 

6.  Multiple  assignment. 

7.  Rewriting  by  the  Laboratory  staff  of  the  A/D  -  D/a  routine 
resulting  in  the  SUBROUTINE  ADDA. 

8.  Rewriting  of  the  digital  processor  control  program 
GATED  to  yield  a  faster  version  GATD1 . 
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9.  The  use  of  equivalencing  scalars  to  vector  elements  to 

reduce  the  indexing  time  within  the  digital  computer. 
10.  Judicious  programming  and  a  dedicated  cleaning  up  process 
of  the  generated  code  for  the  computational  loop. 

The  total  time  for  the  computational  loop  using  the  above  procedures 
was  cut  from  an  original  190  milliseconds  to  55  milliseconds.  This  gives 
an  approximate  solution  update  rate  and  graphic  frame  rate  of  18  samples 
or  frames  per  second. 

The  simulator  was  finally  flown  but  the  dynamic  response  was  not 
quite  as  expected.  The  aircraft  seemed  to  be  unstable  in  the  Dutch  Roll 
mode,  so  a  yaw  damper  was  constructed.  This  helped,  but  at  times  the 
aircraft  was  still  unstable.  Investigation  into  dynamic  error  analysis 
[Ref.  5]  demonstrated  that  the  computational  loop  time  delay  was  causing 
a  degradation  in  the  damping  ratio.  The  lowest  damping  ratio  was  in  the 
Dutch  Roll  mode,  hence  the  appearance  of  the  instability  there.  The 
solution  was  to  predict  ahead  a  value  to  be  used  by  the  program  using 
the  known  derivatives.  When  this  prediction  and  update  was  incorporated, 
all  instabilities  were  corrected. 

Engine  dynamic  response  was  added  in  an  effort  to  make  the  throttle 
response  realistic  [Ref.  6].  Also,  in  an  attempt  to  increase  the  number 
of  cues  available  to  the  pilot,  variable  frequency  noise  was  generated 
in  proportion  to  the  power  applied.  When  weighed  against  the  results 
achieved,  it  was  decided  to  discard  the  noise  generator  in  favor  of  the 
simplicity  of  setup  and  operation. 

The  digital  program  in  its  entirety  can  be  found  in  Appendix  A, 
while  the  analog  diagram  can  be  found  in  Appendix  B.  In  an  effort  to 
make  the  program  more  readable,  Appendix  C  contains  descriptors  of  all 
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the  program  variables  used,  and  Appendix  D  contains  a  program  flow 
diagram.  Sample  output  of  the  program  can  be  found  in  Appendix  G. 
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Simulator 

Aircraft 

.25  cps 

.25  cps 

2.0  sec 

2.3  sec 

37  sec 

38  sec 

VII.  CONCLUSIONS 

Dynamic  validation  of  the  simulator  was  done  by  comparing  the  various 
modes  of  response.  The  tabulation  below  shows  that  in  fact  the  dynamic 
response  of  the  simulator  is  excellent  when  compared  with  the  actual 
aircraft,  in  this  case  the  A-7. 

1.  Dutch  Roll  Natural  Frequency 

2.  Short  Period 

3.  Phygoid  Period 

4.  Roll  Mode  Heavily  damped  in  both  cases 

5.  Spiral  Divergence  Time  Constant  29  sec         25-35  sec 
The  simulator  has  now  met  the  self-imposed  requirements  set  forth 

in  the  Introduction.  Briefly  reviewing  the  points  in  the  Introduction: 

1.  Multi-task  simulator  -  Currently  three  follow-on  studies  are 
being  conducted  using  the  simulator.  These  include  a  spin 
simulation,  investigation  into  stall  dynamics  and  a  study 
involved  with  a  LSO  work  station. 

2.  Ability  to  change  the  aircraft  being  simulated  -  Demonstrations 
have  been  conducted  with  the  "standard"  Navy  runway  and  the 
carrier  runway  dimensions  as  well  as  linearized  and  non- 
linearized  aerodynamic  data.  Time  involved  to  change,  in  each 
case,  was  approximately  one  minute. 

3.  The  use  of  interactive  graphics  for  the  display  -  Very  good 
results  were  achieved  using  the  inside-out  display.  Results 
show  that  there  is  a  learning  curve  as  with  any  simulator  but 
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sufficient  visual  cues  are  provided  to  enable  consistent  landings 
by  experienced  pilots. 

4.  Simplicity  of  setup  and  operation  -  The  number  of  dials  or 

buttons  have  been  kept  to  a  minimum;  as  a  result  the  system  can 
be  up  and  operating  within  five  minutes.  As  for  the  operation, 
the  simulator  has  been  flown  by  all  manner  of  people  with  little 
or  no  explanation  needed. 

Possible  areas  of  follow-on  study  might  include  the  following: 

1.  The  study  of  Heads  Up  Displays  (HUD's),  since  only  the  graphic 
presentation  need  be  changed. 

2.  Autopilot/Gust  Alleviators/Stability  Augmentor  systems  that 
can  be  modeled  on  an  analog  board  can  be  incorporated,  as  approximately 
one-half  the  board  is  left  unused  for  this  purpose. 

3.  Force-stick  analyzation  with  a  realistic  tracking  problem. 

4.  Re-design  of  the  Optical  Landing  System  to  reduce  the  false 
motion  cues  that  currently  exist  in  the  present  system. 
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APPENDIX  A  -  THE  DIGITAL  PROGRAM 

This  appendix  includes  listings  of  the  main  program,  written  in 
FORTRAN,  and  the  SUBROUTINE  ADDA,  the  A/D  -  D/A  subprogram.  The  sub- 
routine is  written  in  METASYMBOL,  the  assembly  language  for  the 
XDS  -  9300  and  was  prepared  gy  the  computer  laboratory  staff. 
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APPENDIX  B  -  THE  ANALOG  PROGRAM 


(R+PQ)    POOO 


Q/SF5 


(P-QR)   P006 
SF5SF6 


Ch 


-PQ    P007 
SF4SF5 


R/SF6 


P/SF4 


■P/SF4 


74 


P/SF4 


MOOO 
N.I. 


Q/SF5 


MOOl 
N.I. 


R/SF6 


PQ/SF4SF5 


-PQ/SF4SF5 


QR/SF5SF6 


M002 
N.I. 


QR/SF5SF6 


R/SF6 


PR 
SF4SF6 


75 


P/SF4 


2 

M015 
D.S. 

x 

P  /SF27        j 

2                * 

-RVSF29 

( 

R/SF6 

1 

2  2 

fr-r 

SF27 


P045 


P011 


-P/SF27 


(P-QR) 
SF5SF6 


-(P-QR) 
SF5SF6 


P012 


-R/SF27 


(R+PQ) 
SF4SF5 


-(R+PQ) 
SF4SF5 


P016 


SX/SF1 


P046 


(R2-P2)/SF27 


(R2-P2)SF27 


76 


P017 


T517 


Sy/SF2 


-S  /SF14 

A021     -3- 


£ 


P020 

O 


-1 


S  /SF3 


P023 


O 


•S  /SF15 


T520 


T422    .  P026 

I I-3/SF6 


T501 


6/SF9 


T504 


4>/SF12 


77 


T503 


e/SFll 


T500 


A045       a/SF8 


T423 


■i/>/SF6 


T502 


^/SFIO 


T420 


V/SF7 


-V/SF1 


78 


-V/SF1 


-V2/(SF1)V^N-PV2/SF20 


O 


P033 


M003  N.I. 


DPDT  SWITCHES 


/wT\ 


-R-PQ/SF4SF5 


QR/SF5SF6  S~*\ 


to  A000 


NO 

v 


-QR/SF5SF6 


to  A000 


o 


2     r2 


IT-PVSF27       Kp\ 


PR/SF4SF6 


■R2+P2/SF27 


0 


/T\ 


to  A001 


-PR/SF4SF6 


NO 
01 


NO 


to  A002 


P-QR/SF5SF6  /^\ 
r    nip     \ 


■P+QR/SF5SF6 


NC 


PQ/5F4SF5  S^\ 
'   NC  ^ 


to  A002 


-PQ/SF4SF5 


02 


0 


NO 


to  A006 


79 


T101 

TM 

-O 

P013 

-1  — 

T510 


TM/SF7 


T105 

Sit 

-o 

P015 

+1- 

T106 


T511 


6it/SF16 


s— O 


+1 


P044 


T512 


6a/SF17 


T104 


T513 


5r/SF18 


-6r/SF18 


80 


T107 


□ 


-1 


P017 


COMPARATORS 


T102 


COO 


00 


-S  /SF13 


■S/SF15 


-1 


Till 


P050 


+1 


01 


; 


C01 


02 


C02 


N422 


C03 


03 


T101 


□ 


P017  SETLINES 

TOOO  ROO  T101         R01 


R02 

T102 

R06 

T103 

1 

UO 

81 


T104  R12 


T105 


R14 


LOO 

• 

LOl 

COMPUTE  R16 


L02 

BINARY  COUNTER 


MC  MCIN  1    kc 

a  n-n 


a  n  6 

CLK  F120  F121  F122 


F123 


dTjtj 


COT  F124  F125  F126  F127 


DF12 


DF12  ENB 


1   usee 


T104 

DF13 

ENB 

DF13 

1  msec 

CLK  N406 


82 


POWER  FOR  LIGHTS 


r\ 


NC 
0 

T112 

+  28  v 

I   no; 

rnc\ 

R13 

+  28  v 

0 

Til  4 

j 

r\ 

R16 

NC 
0 

T113 

+  28  v 

R14 


83 


APPENDIX  C  -  PROGRAM  VARIABLES 

The  following  is  an  alphabetic  listing  of  the  program  variables  used 
in  the  simulator  program. 

A  -  Array  of  angles,  in  order,  alpha,  beta,  psi ,  theta,  and  phi 
A10  -  Initial  conditions  on  A(l),  the  angle  of  attack 
AA  -  Array  of  analog  to  digital  variables 
AD  -  Array  of  analog  to  digital  variables 
ADOTB  -  Normalized  time  derivitive  of  the  angle  of  attack 
ADOTN  -  Array  of  the  negatives  of  the  time  derivitatives  of  the  angles 

in  A 
ADOTO  -  Temporary  storage  for  ADOTN  during  the  predict  and  update  phase 

of  the  program 
ALI  -  Aerodynamic  rolling  moment  about  the  body  axis 
ALPHA  -  Real  argument  based  on  angle  of  attack  and  used  for  coefficient 

lookup 
ALSI  -  Aerodynamic  rolling  moment  about  the  stability  axis 
AMI  -  Aerodynamic  pitching  moment  about  the  body  axis 
ANI  -  Aerodynamic  yawing  moment  about  the  body  axis 
ANSI  -  Aerodynamic  yawing  moment  about  the  stability  axis 
AOA  -  Actual  angle  of  attack  in  degrees 
AOAF  -  Angle  of  attack  break  point  for  fast  cheveron 
AOALDG  -  Optimum  landing  angle  of  attack 

AOALF  -  Angle  of  attack  break  point  for  fast  cheveron  and  donut 
AOALS  -  Angle  of  attack  break  point  for  slow  cheveron 
ARG  -  Argument  used  during  sine  and  cosine  lookup  (type  real) 
ASZ  -  Absolute  value  of  SZ 
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B  -  Wing  span 

B2  -  One  half  the  wing  span  B 

B3  -  Scale  factor  divisor  used  to  divide  out  the  scale  factor  of  the 

homogeneous  transformation  coordinates 
BETA  -  Real  argument  based  on  side  slip  angle  used  in  coefficient 

lookup 
C  -  Array  of  aerodynamic  coefficients  for  a  given  angle  of  attack  and 

sideslip  angle 
CA  -  Array  of  cosines  of  the  angle  array  A 
CB  -  Mean  aerodynamic  chord 
CB2  -  One  half  the  mean  aerodynamic  chord  CB 
CCA(n)  -  Scalar  equivalents  of  the  array  CA 
CHOYC  -  Choice  of  program  options  offered  by  main  program 
COEFA  -  Array  of  coefficients  dependent  on  angle  of  attack  only 
COEFAB  -  Array  of  coefficients  dependent  on  angle  of  attack  and  side 

slip  angle 
C0N1  -  Program  constant  equal  to  the  number  of  degrees  in  a  radian 
C0N2  -  Program  constant  equal  to  .2  divided  by  .75 
C0N3  -  Program  constant  used  in  computing  the  depression  angle  of  the 

horizon 
C0N4  -  Program  constant  equal  to  the  sine  of  30  degrees 
C0N5  -  Program  constant  equal  to  the  cosine  of  30  degrees 
C0N6  -  Program  constant  equal  to  the  ratio  of  C0N3  to  C0N5 
D  -  Constant  equal  to  2.5  times  the  sine  of  the  roll  angle  and  used  in 

calculating  the  horizon 
DA  -  Aileron  deflection 
DEGA  -  Degrees  of  angle  of  attack  for  which  COEFA  and  COEFAB  are  tabulated 
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DEGB  -  Degrees  of  side  slip  for  which  COEFAB  are  tabulated 

DEGMAX  -  Max  angle  of  the  visibility  cone  formed  by  the  square  viewing 

window 
DIT  -  Flying  tail  deflection 
DITO  -  Initial  flying  tail  deflection  setting 
DR  -  Rudder  deflection 
F  -  Focal  length  of  viewing  image 
FXAM  -  Aerodynamic  force  in  the  X-direction  divided  by  the  mass  of  the 

aircraft 
FXSM  -  Force  in  X-direction  of  the  stability  axis  divided  by  the  mass 

of  the  aircraft 
FYAM  -  Aerodynamic  force  in  the  Y-direction  divided  by  the  mass  of  the 

aircraft 
FYSM  -  Force  in  the  Y-direction  of  the  stability  axis  divided  by  the 

mass  of  the  aircraft 
FYWM  -  Force  in  the  Y-direction  of  the  wind  axis  divided  by  the  mass 

of  the  aircraft 
FZAM  -  Aerodynamic  force  in  the  Z-direction  divided  by  the  mass  of  the 

aircraft 
FZWM  -  Force  in  the  Z-direction  of  the  wind  axis  divided  by  the  mass 

of  the  aircraft 
G  -  Acceleration  due  to  gravity 
GEE  -  Scaled  acceleration  due  to  gravity 
GS  -  Glide  slope  in  degrees 
GSO  -  Optimal  glide  slope 

H  -  Graphic  transformation  matrix  (H-Matrix)  or  array 
HI  -  First  component  array  of  the  H-Matrix 
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H2  -  Second  component  array  of  the  H-Matrix 

H3  -  Third  component  array   of  the  H-Matrix 

HHnn  -  Scalar  equivalent  of  the  H-Matrix  or  array 

HINSCT  -  Horizontal  intersection  of  a  line,  used  in  the  software  window 

HMIR  -  Height  of  Lens  below  the  landing  surface 

I  -  Integer  variable  usually  a  counter 

IA  -  Do-loop  counter  based  on  angle  of  attack 

IAGN  -  A  counter  to  check  the  number  of  times  a  point  has  gone  through 

the  software  window 
IALPHA  -  Integer  conversion  of  ALPHA 
IAR6  -  Integer  conversion  of  ARG 
IB  -  Do-loop  counter  based  on  the  side  slip  angle 
IBETA  -  Integer  conversion  of  BETA 
ICHOYC  -  Integer  conversion  of  CHOYC 

IDE  -  Graphic  array  containing  the  dynamic  portion  of  the  display 
IDEV  -  Graphic  device  number  (1  or  2) 

IDIR  -  Graphic  block  directory  for  the  graphics  digital  processor 
IDON  -  l\rray   of  packed  points  for  the  donut  portion  of  the  indexer 
IER  -  Error  flag  returned  by  graphic  subroutines 
IFAST  -  Packed  graphic  array  of  the  fast  cheveron 
IPLS  -  Integer  counter  equal  to  I  plus  one 
ISF7  -  Integer  scale  factor  based  on  thrust  divided  by  mass 
ISLO  -  Packed  graphic  array  of  the  slow  cheveron 

ISQ  -  Array/graphics  block  of  fixed  data,  basically  the  square  or  window 
ITDIR  -  Text  directory  for  the  graphics  digital  processor 
ITEXnn  -  Line  of  text  that  is  sent  from  the  digital  computer  to  the 
graphic  digital  processor 
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IWIRE  -  Wire  caught  on  arrested  landing 

IX  -  Operating  values  of  start/end  points  (X-coordinate)  used  in 

software  window 
1X1  -  Equivalent  of  XSTART,  used  in  software  window 
1X2  -  Equivalent  of  XEND,  used  in  software  window 
IY  -  Operating  values  of  start/end  points  (Y-coordinate) ,  used  in 

software 
IY1  -  Equivalent  of  YSTART,  used  in  software  window 
IY2  -  Equivalent  of  YEND  ,  used  in  software  window 
J  -  Do-loop  counter 

KK  -  Array  of  integers  used  in  outputting  coefficients 
KPOT  -  Array  of  integers  used  in  outputing  analog  pot  settings 
LL  -  Landing  result  flag 

NAD  -  Number  of  analog  to  digital  conversion  variables/trunk  lines 
ND  -  Display  flag  based  on  present  condition  of  the  angle  of  attack 

indexer 
NDA  -  Number  of  digital  to  analog  conversion  variables/trunk  lines 
NULL  -  Octal  value  used  to  null  out  a  line  of  text 
P  -  Angular  velocity  of  the  aircraft  about  X-axis 
Pn  -  Array  of  runway  points 

PDOTN  -  The  negative  of  the  time  derivitive  of  P 
POT  -  Array  of  pot  settings  for  the  analog 
PPnn  -  Scalar  equivalents  of  Pn 
PS  -  P  about  the  stability  axis 
PSB  -  PS  normalized 

Q  -  Angular  velocity  of  the  aircraft  about  the  Y-axis 
QDOTN  -  The  negative  of  the  time  derivitive  of  Q 
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QS  -  Q  about  the  stability  axis 

QSB  -  QS  normalized 

QUE  -  Dynamic  pressure 

R  -  Angular  velocity  about  the  Z-axis 

RAMP  -  X-coordinate  of  the  ramp  of  the  runway 

RDOTN  -  The  negative  of  the  time  derivitative  of  R 

RHO  -  Density  of  air 

RIXX  -  Moment  of  inertia  of  the  aircraft  about  the  X-axis 

RIXZ  -  Product  of  inertia  about  the  X  and  Z-axis 

RIYY  -  Moment  of  inertia  of  the  aircraft  about  the  Z-axis 

RMASS  -  Mass  of  the  aircraft 

RS  -  R  about  the  stability  axis 

RSB  -  RS  normalized 

S  -  Wing  area  of  the  aircraft 

SA  -  Array  of  sines  of  the  angle  array  A 

SDOT  -  Array  of  time  derivitives  (velocities)  of  the  distances  in  all 

three  inertial  directions 
SF1  -  Scale  factor  based  on  maximum  velocity  in  X-direction 
SF2  -  Scale  factor  based  on  maximum  velocity  in  Y-direction 
SF3  -  Scale  factor  based  on  maximum  velocity  in  Z-direction 
SF4  -  Scale  factor  based  on  maximum  angular  velocity  in  X-direction 
SF5  -  Scale  factor  based  on  maximum  angular  velocity  in  Y-direction 
SF6  -  Scale  factor  based  on  the  maximum  angular  velocity  in  the  Z-direction 
SF7  -  Scale  factor  based  on  thrust  divided  by  mass  of  the  aircraft 
SF8  -  Scale  factor  based  on  maximum  angle  of  attack 
SF9  -  Scale  factor  based  on  the  maximum  side  slip  angle 
SF10  -  Scale  factor  based  on  the  maximum  yaw  angle 
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SF11  -  Scale  factor  based  on  the  maximum  pitch  angle 

SF12  -  Scale  factor  based  on  the  maximum  roll  angle 

SF13  -  Scale  factor  based  on  the  maximum  distance  in  the  X-direction 

SF14  -  Scale  factor  based  on  the  maximum  distance  in  the  Y-direction 

SF15  -  Scale  factor  based  on  the  maximum  distance  in  the  Z-direction 

SF16  -  Scale  factor  based  on  the  maximum  flying  tail  deflection  angle 

SF17  -  Scale  factor  based  on  the  maximum  aileron  deflection  angle 

SF18  -  Scale  factor  based  on  the  maximum  rudder  deflection  angle 

SF19  -  Scale  factor  based  on  the  density  of  air 

SF20  -  Scale  factor  based  on  SF1  squared  divided  by  SF19 

SF21  -  Ratio  of  SF6  to  SF4 

SF22  -  Ratio  of  SF1  to  SF4 

SF23  -  Ratio  of  SF1  to  SF5 

SF24  -  Ratio  of  SF1  to  SF6 

SF25  -  Ratio  of  SF4  to  SF5 

SF26  -  Ratio  of  SF5  to  SF6 

SF27  -  SF4  squared 

SF28  -  SF5  squared 

SF29  -  SF6  squared 

SF30  -  Feedback  gain  for  yaw  damper 

SF31  -  Scale  factor  based  on  SF7  divided  by  SF1  and  SF5 

SF32  -  Scale  factor  based  on  SF7  divided  by  SF1  and  SF6 

SF33  -  Ratio  of  SF1  to  SF2 

SF34  -  Ratio  of  SF1  to  SF3 

SF35  -  Scale  factor  based  on  ratio  of  SF5  to  SF8 

SF36  -  Scale  factor  based  on  ratio  of  SF6  to  SF9 

SF37  -  Scale  factor  based  on  ratio  of  SF6  to  SF10 
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SF38  -  Scale  factor  based  on  ratio  of  SF5  to  SF11 

SF39  -  Scale  factor  based  on  ratio  of  SF4  to  SF12 

SF40  -  Ratio  of  SF7  to  SF1 

SFA  -  Array  of  scale  factors,  made  up  to  SF8  through  SF12 

SFACTR  -  Scale  factor  of  the  H-Matrix 

SLOPE  -  Slope  of  a  display  line,  used  in  software  window 

SNKRT  -  Sink  rate  of  the  aircraft  in  ft/sec  at  the  conclusion  of  a  run 

SSAn  -  Scalar  equivalents  of  SA 

SX  -  Actual  distance  of  aircraft  in  X-direction  (inertial) 

SXO  -  Initial  condition  on  SX 

SXN  -  Nromalized  Distance  of  aircraft  in  X-direction  (inertial) 

SY  -  Actual  distance  of  aircraft  in  Y-direction  (inertial) 

SYN  -  Normalized  distance  of  aircraft  in  Y-direction  (inertial) 

SZ  -  Actual  distance  of  aircraft  in  Z-direction  (inertial) 

SZO  -  Initial  condition  on  SZ 

SZN  -  Final  cutoff  value  on  SZ  (height  of  eye) 

T  -  Thrust  of  the  aircraft 

TALPHA  -  Tangent  of  the  yaw  angle 

TAU  -  Time  delay  of  the  prediction  and  update 

TM  -  Thrust  divided  by  mass  of  aircraft 

TPn  -  Temporary  point  number  n 

TRIG  -  Array  of  sines  of  various  angles 

TX  -  X-direction  value  used  to  chop  the  runway 

V  -  Velocity 

VO  -  Initial  condition  on  velocity  V 

VDOT  -  Time  derivitive  of  V 

VDOTO  -  VDOT  of  the  iteration  before 
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VINSCT  -  Vertical  intersection  of  a  line,  used  in  window 

W  -  Weight  of  the  aircraft 

Wn  -  Wire  position  on  runway 

WIND  -  Velocity  of  the  wind  in  knots 

X  -  X  value  used  in  software  window 

XEND  -  Array  of  X-coordinate  ending  points  of  display  lines 

XMIR  -  X  position  of  the  lens 

XSTART  -  Array  of  X-coordinate  starting  points  of  display  lines 

XTEMP  -  Temporary  storage  for  X 

Y  -  Y  value  used  in  software  window 

YEND  -  Array  of  Y-coordinate  ending  points  of  display  lines 

YMB  -  Y-coordinate  of  the  "meatball" 

YMID  -  Y  intercept  of  the  horizon 

YO  -  H-Matrix  offset  in  Y-di recti  on 

YSTART  -  Array  of  Y-coordinate  starting  points  of  display  lines 

YTEMP  -  Temporary  storage  for  Y-coordinate  in  software  window 

Z  -  Temporary  storage  array  used  in  multiplying  HI,  H2,  and  H2  to  yield 

the  H-Matrix 
ZO  -  H-Matrix  offset  in  Z-direction 
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APPENDIX  D  -  DIGITAL  PROGRAM  FLOW  DIAGRAM 
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APPENDIX  E  -  OPERATING  MANUAL 

The  program  exists  on  tape  in  two  forms,  one  is  in  a  Source  deck 
format  (SI)  and  must  be  compiled  prior  to  execution.  This  takes  approxi- 
mately 15  minutes.  This  tape  is  stored  in  the  tape  locker  and  is  titled 
CALS  NEW  VERSION.  The  other  version  is  a  core  dump  (SAVE)  version.  Just 
prior  to  the  execution  of  the  program,  after  the  compilation  and  loading 
phases  are  complete,  a  copy  of  core  is  dumped  on  tape.  Execution  of 
this  version  (RERUN)  is  approximately  10  seconds.  Again  this  is  stored 
in  the  tape  locker  under  the  title  of  CALS  CORE  DUMP. 

All  the  required  cards  for  the  program  are  filed  in  the  file  under 
CARRIER  APPROACH  LANDING  SIMULATOR.  The  cards  filed  include  the  master 
version  of  the  program,  several  data  decks  and  the  control  cards  for 
running  any  of  the  different  programs.  The  control  cards  are  color 
coded  with  explanations  of  their  use  written  on  each  set. 

The  instructions  for  operating  the  digital  computer  and  graphics 
display  are   already  published  and  can  be  obtained  from  the  Laboratory 
Office. 

The  following  is  a  sequence  of  events  and  instructions  for  operating 
the  simulator  assuming  all  equipment  is  turned  on. 

1.  Load  the  version  of  tape  selected  on  the  tape  drives.  Select 
the  proper  mount  (2  for  SI,  3  for  RERUN). 

2.  Load  and  ready  the  card  reader  with  control  cards  followed 
by  the  DATA. 

3.  Ready  the  Line  Printer. 

4.  On  the  Digital  Main  Frame,  Punch  IDLE,  RESET,  RUN,  CARDS. 
Cards  and  tape  should  start  being  read  at  this  time. 
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5.  Load  the  Analog  Program  board  and  the  Analog  Logic  board 
into  the  computer  (Boards  number  7). 

6.  On  the  Analog  Keyboard  punch  KEYBOARD,  LOCAL,  POTSET.  In  the 
far  right  control  box,  only  the  following  can  be  lit:  IDCX1  or 
IDCX.l  and  REAL  TIME.  All  other  lights  must  be  punched  out. 
When  complete  punch  DIGITAL  CMPTR. 

7.  On  the  Analog  the  following  options  exist: 
DSO  -  UP  -  no  printed  output 

CENTER  -  printed  output 
DS1  -  UP  -  distance  integrators  disabled 

CENTER  -  distance  integrators  enabled 
DS2  -  UP  -  yaw  damper  ON 

CENTER  -  yaw  damper  OFF 
DOWN  positions  on  these  switches  are  momentary  contact  positions, 

8.  On  the  Analog: 

Turn  the  inner  dials  of  DF12  and  DF13  clockwise  until  resistance 
is  felt.  Then  set  the  outer  dial  of  DF12  to  .1  and  the  outer 
dial  of  DF13  to  10. 

9.  Under  Master  Clock  (Analog)  punch  RUN.  The  run  light  on  the 
light  panel  should  be  on  and  the  DF13  light  should  be  flashing. 

10.  At  the  Graphics  terminal  teletype  to  be  used  (usually  number  2), 
load  the  control ing  program  by  typing 

RESET  ("GATD1",  120)  ! 
If  the  above  is  answered  with 

FILE  NOT  FOUND 
type    RESET  ("GATD1",  20)! 

Then  execute  the  program  by  typing  GATED!    At  this  point 
nothing  should  be  on  the  display  screen. 
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11.  By  this  time  the  Digital  computer  should  have  compiled  and 
loaded  the  main  program  and  the  following  will  be  typed  out  on 
the  Digital  consol  teletype: 

INPUT  AGT  NUMBER 
The  input  light  will  then  come  on.  To  respond,  type  the  AGT 
number  (1  or  2)  of  the  unit  being  used,  followed  by  a  carriage 
return. 

12.  The  Title  should  appear  on  the  selected  AGT. 

13.  The  Data  will  be  read  in  from  the  card  reader. 

14.  If  output  is  selected,  the  output  will  be  printed  on  the  line 
printer. 

15.  The  Analog  Pots  will  be  set  by  the  Digital  computer.  The 
Address  and  Ratiometer  readings  should  be  changing. 

16.  Set  the  cockpit  in  front  of  the  AGT  (will  only  reach  AGT-2). 
The  cockpit  is  stored  behind  the  Analog  computer. 

17.  The  Instructions  should  flash  on  the  screen. 

18.  Punch  the  button  on  the  Throttle  plate  and  the  display  will 
appear. 

19.  Punch  the  button  on  the  control  stick  to  fly.   If  during  a  run 
you  wish  to  abort  that  run,  punch  the  button  on  the  Throttle 
plate. 

20.  At  the  completion  of  the  run  the  Landing  Results  will  be 
displayed,  followed  by  a  short  delay  in  which  to  read  them. 

21.  At  this  time  two  options  are  offered.  To  fly  again  punch  the 
button  on  the  Control  stick  and  the  program  jumps  back  to  the 
Instructions  (#17). 
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22.  To  receive  the  expanded  program  options,  punch  the  Throttle 
button.  Follow  the  displayed  instructions  and  type  your 
selection  on  the  AGT  teletype  unit.  This  will  result  in  one 
of  three  things  happening. 

1.  Logical  STOP  to  the  program. 

2.  Jump  to  the  Instructions  (#17). 

3.  Jump  to  the  selection  of  the  AGR  number  (#11). 

23.  At  the  completion  of  the  computer  time.  Return  all  used 
equipment  (cards,  tape,  boards,  cockpit)  to  their  storage 
locations. 
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APPENDIX  F  -  PREPARATION  OF  DATA  DECK 

The  following  is  a  guide  to  be  used  in  preparing  the  data  deck  to 
be  used  with  the  simulator. 

Cards  1  -  6  -  The  first  six  data  cards  consist  of  the  actual 
coordinates  of  the  six  runway  points.  The  first  ten  columns  of  the  card 
are  the  x  coordinate  in  feet,  the  second  ten  columns  are  the  y  coordin- 
ate in  feet,  the  third  ten  columns  are  the  z  coordinate  in  feet  and  the 
next  ten  columns  are  the  H-Matrix  scale  factor  which  in  this  simulator 
is  1.0.  The  remaining  columns  are  left  blank.  The  order  of  the  runway 
points  is  as  follows:  far-centerline,  far-left,  near-left,  near-center- 
line,  near-right,  far-right. 

Card  7  -  The  first  ten  columns  consist  of  the  optimum  glide  slope 
(GSO)  in  degrees.  The  second  ten  columns  contain  the  x  coordinate  of 
the  location  of  the  Fresnel  Lens  (XMIR)  in  feet.  The  third  ten  columns 
contain  the  velocity  of  the  headwind  (WIND)  in  knots. 

Card  8-10  -  These  cards  contain  the  inputed  scale  factors  for  the 
scaling  of  the  analog  solution.  Each  scale  factor  has  a  field  of  ten 
columns  and  are  ordered  as  follows: 

Card  8  -  SF1,  SF2,  SF3,  SF4,  SF5,  SF6,  SF8,  SF9 
Card  9  -  SF10,  SF11,  SF12,  SF13,  SF14,  SF15,  SF16,  SF17 
Card  10  -'  SF18,  SF19,  SF30,  the  remaining  columns  are  blank 
Card  11  -  This  card  contains  the  initial  conditions  of  the  problem 
in  ten  column  increments  as  follows:  initial  x  coordinate  (SXO)  in  feet, 
initial  altitude  (SZO)  in  feet,  initial  velcoity  (VO)  in  ft/sec,  initial 
angle  of  attack  (A10)  in  degrees,  initial  elevator  angle  (DITO)  in 
degrees  and  the  density  of  the  air  (RHO)  in  slug/ft  .  The  remaining 
portion  of  the  card  is  left  blank. 
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Card  12  -  This  card  contains  a  portion  of  the  aircraft  constants 
again  in  ten  column  intervals.  They  are  as  follows:  weight  (W)  of  the 

aircraft  in  lb,  wing  span  (B)  in  ft,  wing  chord  (CB)  in  ft,  wing  area 

2  2 

(S)  in  ft  ,  moment  of  inertia  about  the  x  axis  (RIXX)  in  slug-ft  , 

moment  of  inertial  about  the  y  axis  (RIYY)  in  slug-ft  ,  moment  of  inertia 

p 
about  the  z  axis  (RIZZ)  in  slug-ft  ,  and  the  cross  product  of  inertia 

(RIXZ)  in  slug-ft2. 

Card  13  -  This  card  contains  some  more  aircraft  constants  again  in 

ten  column  intervals.  They  are  as  follows:  maximum  thrust  (T)  in  lbs, 

optimum/landing  angle  of  attack  (AOALDG)  in  degrees,  height  of  eye  at 

touchdown  (SZF)  in  feet.  The  remaining  columns  are  left  blank. 

Cards  14  -  85  (Alternating)  -  These  cards  contain  the  aerodynamic 

coefficients  dependent  on  both  alpha  and  beta.  The  first  card  in  the 

series  is  as  follows:  in  columns  1-2  the  coefficient  number,  columns 

3-4  blank,  column  5  the  index  on  beta  (IB),  colums  6-10  blank.  On  the 

remaining  portion  of  the  card  in  ten  column  intervals  are  the  coefficients 

indexed  on  alpha  (IA).  The  second  card  in  the  series  also  contains  the 

remaining  coefficients  again  in  ten  column  intervals  and  indexed  on  alpha. 

Alpha  is  indexed  15  times  for  each  beta  index  while  beta  is  indexed  9 

times  for  each  of  the  four  coefficients.  The  range  on  alpha  is  -25 

degrees  to  +  45  degrees  with  5  degree  intervals  while  the  range  on  beta 

is  -40  degrees  to  +40  degrees  with  10  degree  intervals. 

The  coefficients  are  as  follows: 

Number 

1 

2 

3 

4 
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Coefficient 

ci 

cm 

Cn 

cY 

Cards  86  -  125  (Alternating)  -  These  cards  contain  the  aerodynamic 
coefficients  dependent  only  on  alpha.  The  first  card  in  the  series  is 
as  follows:  columns  1-2  the  coefficient  number,  columns  3-10  blank.  The 
remaining  portion  of  the  card  in  ten  column  intervals  are  the  aerodynamic 
coefficients  indexed  on  alpha.  The  second  card  of  the  series  contains 
the  remaining  coefficients  again  in  ten  column  intervals  and  indexed  on 
alpha.  Alpha  is  indexed  15  times  for  each  of  the  coefficients.  A  list 
of  the  remining  coefficients  and  their  numbers  is  given  below. 

Number  Coefficient 

5  CD 

6  CL 

7  CL  . 

L6it 

8  Cm 

m6it 

9  CY 

6r 

10  C, 

6r 

11  Cn 

V 

12  CY 

6a 

13  C, 

6a 

14  C 

n6a 

15  CY 

YP 

16  C, 

'P 

17  Cn 

18  C 


YR 


19  C] 
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R 


Number 

Coefficient 

20 

\ 

21 

\ 

22 

C 
mQ 

23 

C 

a 

24 

m 
a 
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APPENDIX  G  -  SAMPLE  OUTPUT 

The  following  is  a  sample  output  of  the  digital  program  and  contains 
the  following  sections: 

1 .  Runway  Constants 

2.  Scale  Factors  and  Aircraft  Constants 

3.  Original  Aircraft  Aerodynamic  Coefficients 

4.  Usable  Aircraft  Aerodynamic  Coefficients 

5.  Pot  Settings  for  the  Analog  Computer 

6.  Output  of  the  Amplifiers  and  D/A  Trunks 
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